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(T) ' We consider a gravastar model made of anisotropic dark energy with an infinitely thin spher- 
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CN ■ ical shell of a perfect fluid with the equation of state p = (1 — j)o~ with an external de Sitter- 

. Schwarzschild region. It is found that in some cases the models represent the "bounded excursion" 
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first time in the literature that a naked singularity emerges from a gravastar model. 
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I. INTRODUCTION 



Nowadays, although we know that gravastars are not an alternative to black holes, these 
theoretical objects cannot be ignored and, in fact, they have received a considerable attention 

The pioneer model of gravastar was proposed by Mazur and Mottola (MM) |4| , consisting 
of five layers: an internal core < R < R± , described by the de Sitter universe, an 
intermediate thin layer of stiff fluid R± < R < R2 , an external region R > R2, described 
by the Schwarzschild solution, and two infinitely thin shells, appearing, respectively, on the 
hypersurfaces R = R\ and R = R 2 , where R denotes the radius of the star. The intermediate 
layer is constructed in such way that R\ is inner than the de Sitter horizon, while R 2 is outer 
than the Schwarzschild horizon, eliminating the apparent horizon. Configurations with a de 
Sitter interior have long history which we can find, for example, in the work of Dymnikova 
and Galaktionov After this work, Visser and Wiltshire j(3] pointed out that there are 
two different types of stable gravastars which are stable gravastars and "bounded excursion" 
gravastars. In the spherically symmetric case, the motion of the surface of the gravastar can 



be written in the form 



6], 



l -R 2 + V(R) = 0, (1) 



where R = dR/dr, with r being the proper time of the surface. Depending on the properties 
of the potential V(R), we can have two kinds of gravastars, which are (i) Stable gravastars 
and (ii)" Bounded excursion" gravastars. For the former one there must exist a radius Rq 
such that V (Rq) = 0, V (Rq) = 0, V" (Rq) > 0, where a prime denotes the ordinary 
differentiation with respect to the indicated argument. If and only if there exists such a 
radius Rq for which the above conditions are satisfied, the model is said to be stable. 

On the other hand, for the second one, there exist two radii Ri and R2 such that V (Ri) = 
0, V (Rx) <0,V (R 2 ) = 0, V (R 2 ) > 0, with V(R) < for a G (R u R 2 ), where R 2 > R x . 

Many authors have studied the stabilitv of the gravastar models for several equations of 
state, among them we can mention f| - 9]. Some generalizations of these models can be 
found in the literature -0. 

The study of gravastar, in general, has considered these objects embedded in a 
Schwarzschild spacetime (except in the references [7( and [HI]). However, taking the cosmo- 



logical point of view that the universe must be fulfilled by a considerable amount of dark 



2 



energy, it is very important to investigate its influence in the gravastar stability and in its 
possible dynamic evolution. In a first step, we have considered the de Sitter-Schwarzschild 
exterior spacetime, in order to introduce a positive cosmological constant, which has been 
suggested SIS db dark energy candidate. 



In this paper, we generalize our previous works ll|-[13| to the case where the equation 



of state of the infinitely thin shell is still given by p = (1—7)0" with 7 being a constant, 



the interior consists of an anisotropic dark energy fluid (similarly used in 13]), while the 



exterior is now the de Sitter-Schwarzschild space (similarly used in [ll|]). In the previous 



work 13| we showed that anisotropic dark energy could collapse forming a black hole or a 
gravastar. We also studied the effect of anisotropy of the interior fluid on the formation of 
gravastars and it was concluded that the sign of the difference between the pressures (radial 
and tangential) affects the conditions of the formation of the gravastar and black holes when 
the interior fluid of prototype gravastars is anisotropic. This is confirmed in this present 
work. 

Here we shall first construct three-layer dynamical models, and then show both types of 
gravastars and black holes exist for various situations. In addition, this model shows that a 
naked singularity can be the final stage of a gravastar. The rest of the paper is organized 
as follows: In Sec. II we present the metrics of the interior and exterior spacetimes, and 
write down the motion of the thin shell in the form of equation ([1]). In Sec. Ill we show the 
definitions of dark and phantom energy, for the interior fluid and for the shell. In Sec. IV 
we discuss the different structures that are formed from standard or dark energy anisotropic 
fluid. Finally, in Sec. V we present our conclusions. 

II. DYNAMICAL THREE-LAYER PROTOTYPE GRAVASTARS 

The interior fluid is made of an anisotropic dark energy fluid with a metric given by the 
first Lobo's model [3] 

ds 2 = -f\dt 2 + f 2 dr 2 + r 2 dQ 2 , (2) 
where dVL 2 = d6 2 + sin 2 (0)# 2 , and 

,2\ 



h = (l-2ar 2 )-— , 
h = 1^2^' (3) 
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where u is a constant, and its physical meaning can be seen from the following equation 
(jSJ)- Since the mass is given by m(r) = 47rpor 3 /3 and a = 4irpo/3 then we have that a > 0, 
where po is the homogeneous energy density. Note that there is a horizon at = 1/ y2a, 
thus the radial coordinate must obey r < r^. The corresponding energy density p, radial 
and tangential pressures p r and pt are given, respectively, by 



P 

Pr 
Pt 



p = constant, 



wp 



1 



4tt(1 +w)(l + 3w)p r 2 



6 u(l 



(4) 
(5) 



when a; = — 1 and u = -1/3 we obtain an interior isotropic pressure fluid. 
The exterior spacetime is given by the de Sitter-Schwarzschild metric 



(6) 



ds 2 e = -fdv 2 + /"Mr 2 + rW, 

where f = 1 — ^ — (r/L e ) 2 and L e = a/3/A 6 . The metric of the hypersurface on the shell 
is given by 

dsl = -dr 2 + R 2 (r)dn 2 , (7) 

where r is the proper time. 

Since ds 2 = dsl = ds 2 ^, we find that = = R, and 



fit 2 - / 2 r| = 1, 
r 2 

'*-7 = 1 ' 



(8) 
(9) 



where the dot denotes the ordinary differentiation with respect to the proper time. On the 
other hand, the interior and exterior normal vectors to the thin shell are given by 



n. 



;-r s ,i,0,0), 
:-r s ,i;,0,0). 



(10) 



Then, the interior and exterior extrinsic curvatures, using the junction condition = 
R and equations ([8]), (|9]) calculated at the shell radius, are given by 



& = _(i _ 2aj R2 ) -(3a, + i)/ 2 f f 6(1 _ 2a R 2 )^ +1 ^ 2 R 2 uj + QaR 2 i 2 u + 2aR 2 i 2 - 3t 2 uj - i 2 



:aRi - (1 - 2a J R 2 ) (3w+1)/2 (-l + 2aR 2 )(Rt - Ri)} (-1 + 2aR 2 y 1 . i 



(11) 



K\ e =i(l- 2aR 2 )R, 
K^ = K l ee sin 2 (6), 

K e TT = v[{2L 2 e mv + L 2 RR - L 2 e Rv + R 3 v)(2L 2 e mv - L\RR - L 2 Rb + R 3 v) - 
2L 4 e R 2 R 2 }((2m - R)L\ + R 3 )-\L 2 e m - R 3 )L^R~ 3 + Rv - Rv 

K e ee = -v((2m-R)L 2 e + R 3 )L; 2 
K e u = K e 0e sm 2 (9). 



Since |20| 



[K e 



^00 ^0 



-M, 



where M is the mass of the shell, we find that 



M = vR 



1 



2m 



R 



R \ L. 

Then, substituting equations dSJ) and (EJ) into (|T8|) we get 



+ t(l-2aR 2 )R. 



(12) 
(13) 

(14) 

(15) 
(16) 

(17) 
(18) 



M 



-R 



2m 



R_ 



+ R 2 



1/2 (l - 2aR 2 + R 2 ) 
+ R K 1 



1/2 



[1 - 2aR 2 )~(^+ 2 )/ 2 



(19) 



In order to keep the ideas of MM as much as possible, we consider the thin shell as consisting 
of a fluid with the equation of state, p = (1 — 7)0", where a and p denote, respectively, the 
surface energy density and pressure of the shell and 7 is a constant. Then, the equation of 



motion of the shell is given by 



2Q| 



M + SttRRp = 4nR 2 {T af3 u a nP} = AtiR 2 (T^u^ e - 3^<nf) , (20) 
where u a is the four-velocity. 

P a/3 U 'e ,v e ~ - L a(3 u i 'H 



Let us first calculate T^w^nf — T^ufn 1 ^. We can get the energy-momentum tensor from 



reference 



191 ] given by 



T aP = (P + Pt)K u p + PtQaP + (Pr - Pt)XaXf), 



(21) 



where u l a is the timelike four- velocity vector (u l a Ui a = —1), Xa is the unit spacelike vector 
in the radial direction (xaX a = 1)> n a * s the spacelike normal vector n l a rii a = jj-^ (fx > 
and f% > 0), in the region filled with the anisotropic fluid. Thus, we have 



<x a = 0, 



(22) 



and 

= 0. (23) 
Thus, since = 0, we can easily proof that T^u^n^ — T^u^n^ = 0. Thus, 

M + 87^(1- 7)0- = 0. (24) 

Recall that a = M/(4tcR 2 ), we find that equation ( I24p has the solution 

M = kR 2 ^~ 1 \ (25) 

where k is an integration constant. Substituting equation (|25|) into equation (TT9]) . and 
rescaling m, L e a and R as, 

1 

m — > mk 2 t- 3 , 
2 

a — > ak 2 ~/- 3 , 
L e L e k^, 

R Rk~^, (26) 
we find that it can be written in the form of equation ([1]) with a replaced by R, and 



V{R, m, L e , u, a, 7) 
1 



2R 2 ^ 2 bL e (wzA {R 2 ^ 2 bL e 



2R 2 L 2 e (b 2 - : 

^2mL\R + R 4 - 2R 4 L 2 e a + L\R^ - 2b 2 mL 2 e R - b 2 R 4 + 2R i L 2 e ab 2 ' 
+R 2 L 2 e - R 2 L\b 2 + 2R i L 2 e ab 2 - 2mL 2 e R - R 4 - L 2 e R^\ (27) 



where 



b= (l-2ai? 2 )- (1+3w/2) . (28) 

Clearly, for any given constants m, w, a and 7, equation (|2"T|) uniquely determines the collapse 
of the prototype gravastar. Depending on the initial value Rq, the collapse can form either a 
black hole, a gravastar, a de Sitter, or a spacetime filled with anisotropic dark energy fluid. 
In the last case, the thin shell first collapses to a finite non-zero minimal radius and then 
expands to infinity. 



The exterior horizons are given by 23 1 



2m ( 7r + ib\ 
r bh = —cos — ^ , (29) 



3y 



rc = 7¥ cos V^J' m 

where y = (m/L e ) 2 , ip = arccos (3y/3y), r bh denotes the black hole horizon and r c denotes 
the cosmological horizon. 

To guarantee that initially the spacetime does not have any kind of horizons, cosmological 
or event, we must restrict Ro to the range, 

r bh < Ro <r h or r c , (31) 

where Rq is the initial collapse radius. When m = = a, the thin shell disappears, and the 
whole spacetime is Minkowski. So, in the following we shall not consider this case. 

Since the potential ( 127|) is so complicated, it is too difficult to study it analytically. 
Instead, in the following we shall study it numerically. 



III. CLASSIFICATIONS OF MATTER, DARK ENERGY, AND PHANTOM EN- 
ERGY FOR ANISOTROPIC FLUIDS 

Recently |22j, the classification of matter, dark and phantom energy for an anisotropic 
fluid in terms of the energy conditions was studied since the pressure components may play 
very important roles and can have quite different contributions. In this paper, we will use 
this classification to study the collapse of the dynamical prototype gravastars, constructed 
in the last section. The denomination used for each kind of fluid is given in Table 1. 

In order to consider the equations ([2]) and ([5]) for describing dark energy stars we must 
analyze carefully the ranges of the parameter u that in fact furnish the expected fluids. It 
can be shown that the condition p + p r > is violated for u> < — 1 and fulfilled for ui > —1, 
for any values of R and b. The conditions p + Pt > and p + p r + 2p t > are satisfied 
for uj < — 1 and —1/3 < u < 0, for any values of R and b. For the other intervals of u 
the energy conditions depend on very complicated relations of R and b [22]. This provides 
an explicit example, in which the definition of dark energy must be dealed with great care. 



Another case was provided in a previous work [22 ]. 

In order to fulfill the energy condition a + 2p > of the shell and assuming that p = 
(1 — 7)0" we must have 7 < 3/2. On the other hand, in order to satisfy the condition 
cr + p > 0, we obtain 7 < 2. Hereinafter, we will use only some particular values of the 
parameter 7 which are analyzed in this work. See Table II. 



TABLE I: This table summarizes the classification of the interior matter field, based on the energy 
conditions 211 ]. where we assume that p > 0. 



Matter Condition 1 Condition 2 Condition 3 
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TABLE II: This table summarizes the classification of matter on the thin shell, based on the 

™ s . COndittos fl. Th e ,a, t C0tam „ ^ th e P a rticular ^ „ f the p_ , w here 
we assume that p > 0. 

Matter Condition 1 Condition 2 7 

Standard Matter a + 2p>0 a + p>0 -lorO 

Dark Energy a + 2p < a+p>0 7/4 

Repulsive Phantom Energy cr + 2p<0 <7+p<0 3 



In the next sections we will discuss the different possibilities for the type of system that 
can be formed from the study of the potential V{R 1 m 1 L ei u 1 a 1 r ))\ black hole of standard 
matter or dark energy, stable or "bounded excursion" gravastar and even naked singularity. 



IV. STRUCTURES FORMED 



Here we can find many types of systems, depending on the combination of the constitution 
matter of the shell and core. Among them, there are formation of black holes, stable and 
"bounded OTi on" gra v as ta r8 , as it has already snow, in on r prev.ons w„ rk s fl-Q. and 
even a naked singularity constituted exclusively of dark energy. All of them are listed in the 
table III. 



S 



As can be seen in the figure [Q depending on the value of the cosmo logical constant, 
we can see that V(R) = now can have one, two or three real roots. Then, we have, 
say, Ri, where R i+ i > Ri. For L e = L\ (corresponding to A = 0) If we choose _R > R 3 
none structure is allowed in this region because the potential is greater than the zero. 
However, if we choose R2 < Ro < R3, the collapse will bounce back and forth between 
R = Ri and R = i?2- Such a possibility is better shown in the figure [9j This is exactly 
;he so-called "bounded excursion" model mentioned in [6|, and studied in some details in 
8j- 12|. Of course, in a realistic situation, the star will emit both gravitational waves and 
particles, and the potential will be self- adjusted to produce a minimum at R = Rstntir. where 
V (R = Rstatic) = = V (R = Rgtatic) whereby a gravastar is finally formed 10|. For 
Ro < R\ a black hole is formed in the end of the collapse of the shell. 

The scenario above can significantly be changed if we consider A > 0. In this case for 
L e > L c , we also have bounded excursion gravastars if R 2 < Ro < R-s- However, for R < R± 
the final structure can be now a black hole or a naked singularity since the presence of the 
cosmological constant above a certain limit (L*) eliminates the event horizon (its radius 
becomes negative), as can be seen in the table IV. See figures HJ [5] and El This is the first 
evidence of a naked singularity formation from a gravastar model. Moreover for L e = L c , 
then i?2 = -R3, a stable gravastar is formed if Ro = R2, while for L e < L c there is only one 
real root. Note that for any value of L e > L*, a naked singularity is formed for small initial 
radius of the shell. 

Another two very interesting structures that arise from this model are a black hole, 
represented by figures dOj [15] and HEl and a stable gravastar, represented by figure [HI 
all of them made of dark energy. It is remarkable that the stable gravastar is constituted 
only by phantom energy. It means that a system constituted only by dark energy is able to 
achieve an equilibrium state or even to collapse. Thus, solving equation ( TT9|) for R(t) we 
can integrate R(t) and obtain R(t), which are shown in the figures El El [TT1 and [T2l for the 
case F. 



V. CONCLUSIONS 

In this paper, we have studied the problem of the stability of gravastars by constructing 
dynamical three-layer models of VW [6] , which consists of an internal anisotropic dark energy 



9 



TABLE III: This table summarizes all possible kind of energy of the interior fluid and of the shell 
and compares the formed structures in the two gravastar models (A = 0, A > 0). The letters 
SG, UG, BEG, BH, NS and N denote stable gravastar, unstable gravastar, bounded excursion 
gravastar, black hole, naked singularity and none, respectively. 



Case Interior Energy 



Shell Energy 



Figures Structures (A = 0) Structures (A > 0) 



A Standard Standard BH BH 

B Standard Dark BH BH 

C Standard Repulsive Phantom UG/BH BH 

D Dark Standard SG N 

E Dark Dark M UG/BH BH 

F Dark Repulsive Phantom HEgEEEEl BH/BEG BEG/SG/NS 

[T0TT31 BH/UG BH 

G Repulsive Phantom Standard BH BH 

H Repulsive Phantom Dark QU BH BH 

I Repulsive Phantom Repulsive Phantom [J3] SG N 

J Attractive Phantom Standard N N 

K Attractive Phantom Dark N N 

L Attractive Phantom Repulsive Phantom N N 



fluid, a dynamical infinitely thin shell of perfect fluid with the equation of state p = (1—7)0", 
and an external de Sitter-Schwarzschild spacetime. 

We have shown explicitly that the final output can be a black hole, a "bounded excursion" 
stable gravastar depending on the total mass m of the system, the cosmological constant 
L e , the parameter u, the constant a, the parameter 7 and the initial position Rq of the 
dynamical shell. All these possibilities have non-zero measurements in the phase space of 
m, L e , a, u, 7 and Rq. All the results can be summarized in Table III. 

An interesting result that we can deduce from Table III is that we can have black hole 
and stable gravastar formation even with an interior and a shell constituted of dark and 
repulsive dark energy (cases E, F, H and I). 

We also would like to point out the significant influence of the presence of the exterior 
cosmological constant to formation of this kind of structure. Note that in the case I, rep- 
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TABLE IV: Radii rt,h and r c calculated using equation (|29p as a function of L e for L e < L c and 



for L e > L c where m = 10~ 9 and L c = 61.49177124. 
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FIG. 1: The potential V{R) and the energy conditions EC1= p + p r + 2pt, EC2= p + p r and 
EC3= p+pt, for 7 = 3, u = -2/3, a = 0.01 and m c = 0.1054688609 x 10" 8 . The horizons are: r h = 
7.071067814; r bh = -0.2368329842 x 10~ 7 , r c = 99.99999998 (L e = 100); r bh = -0.1657830889 x 
10- 7 , r c = 70 (L e = 70); r bh = -0.1456327969 x 10~ 7 , r c = 61.49177126 (L e = 61.49177124). Note 
that we have a gravastar enclosing a naked singularity. Then, for small R, the shell can collapse 
to form a naked singularity. Case F 

resented by figure dU we have a stable gravastar (for A = 0) or no formed structure (for 
A > 0), as well as in the case D, where we have a traditional stable gravastar only for A = 0. 
Still more interesting is the case F, represented by figure [TJ where for small radius of the 
shell we have formation of a black hole (for A = 0) or a naked singularity (for A > 0). 
This is the first time in the literature that a naked singularity emerges from a gravastar 
model. Besides, the figures El El [11] and [12] give us examples of the dynamical evolution of 
a gravastar to a naked singularity, a "bounded excursion" gravastar, a black hole and an 
unstable gravastar, respectively. 
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FIG. 2: The zoom of the potential V(R) (figured]) and the energy conditions EC1= p + p r + 2pt, 
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(L e = 61.49177124). Note that we have a gravastar enclosing a naked singularity. Then, for small 
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FIG. 6: The zoom of the potential V(R) (figured]) and the energy conditions EC1= p + p r + 2pt, 
EC2= p + pr and EC3= p + p t , for 7 = 3, u = -2/3, a = 0.01 and m c = 0.1054688609 x 10~ 8 . 
The horizons are: r h = 7.071067814; r bh = -0.1184164921 x 10~ 7 , r c = 50.0 (L e = 50); r bh = 
-0.9473319370 x 10" 7 , r c = 40 (L e = 40); r bh = -0.7104989528 x 10~ 8 , r c = 30.0 (L e = 30). Note 
that we have a gravastar enclosing a naked singularity. Then, for small R, the shell can collapse 
to form a naked singularity. Case F 
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FIG. 7: The zoom of the potential V(R) (figured]) and the energy conditions EC1= p + p r + 2pt, 
EC2= p + p r and EC3= p + pt, for 7 = 3, u = -2/3, a = 0.01 and m c = 0.1054688609 x 10~ 8 . The 
horizons are: r h = 7.071067814; r bh = 0.9178675538 x 10~ 8 , r c = 10.0 {L e = 10). Note that we 
have a gravastar enclosing a black hole. Then, for small R, the shell can collapse to form a black 
hole. Case F 




FIG. 8: These figures represent the dynamical evolution of a "bounded excursion" gravastar with 
the potential given by the figure [H assuming Rq = R(0) = 0.007. 
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FIG. 9: Zoom of the potential V(R) of figure Q] and the energy conditions EC1= p + p r + 2pt, 
EC2= p + pr and EC3= p + p t , for 7 = 3, u = -2/3, a = 0.01 and m c = 0.1054688609 x 10~ 8 . 
The horizons are: r h = 7.071067814; r bh = -0.2368329842 x 10~ 7 , r c = 99.99999998 (L e = 100); 
r bh = -0.1657830889 x 10" 7 , r c = 70 (L e = 70); r bh = -0.1456327969 x 10" 7 , r c = 61.49177126 
(L e = 61.49177124). Case F 
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FIG. 10: The potential V(R) and the energy conditions EC1= p + p r + 2pt, EC2= p + p r and 
EC3= p + p t , for 7 = 3, w = -2/3, a = 0.01 and m c = 0.1054688609 x 10" 8 . The horizons 
are: r h = 7.071067814; r bh = 0.1891503488 x 10~ 8 , r c = 0.9126375596 (L e = 0.9126375605); 
r bh = 0.1761682879 x 10~ 8 , r c = 0.8499999992 (L e = 0.85). Case F 




FIG. 11: These figures show the dynamical evolution of the collapse of a gravastar with the potential 
given by the figure [T0| assuming Rq = R(0) = 0.9, forming at the end of the evolution a black hole. 
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0.4 




FIG. 12: These figures show the dynamical evolution of an unstable gravastar with the potential 
given by the figure HH1 assuming Rq = R(0) = 0.9. 




FIG. 13: The potential V{R) and the energy conditions EC1= p + p r + 2pt, EC2= p + p r and 
EC3= p + pt, for 7 = 3, uj = —2/3, a = 0.01 and m c = 0.5104531527. The horizons are: = 
7.071067814; r bh = 1.022092711, r c = 29.47589240 (L e = 30); r bh = 1.022617334, r c = 24.47300022 
(L e = 25); r bh = 1.023587410, r c = 19.46855167 (L e = 20). Case F 
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FIG. 14: The potential V(R) and the energy conditions EC1= p+p r + 2pt, EC2= p+p r and EC3= 
p + pt, for 7 = 3, u = -3, a = 0.1 and m c = 0.0479271625. The horizons are: r h = 2.236067977; 
r bh = 0.09588959554, r c = 4.951365546 (L e = 5); r bh = 0.09595248530, r c = 2.950872678 (L e = 3). 
Case I 
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FIG. 15: The potential V(R) and the energy conditions EC1= p+p r + 2p t , EC2= p+p r and EC3= 
p+pt, for 7 = 7/4, u = -2/3, b = 0.01 and m c = 0.79402438. The horizons are: r h = 7.071067814; 
r bh = 1.588449580, r c = 99.19631286 (L e = 100); r bh = 1.589655604, r c = 49.18621608 (L e = 50); 
r bh = 1.592536495, r c = 29.17201284 (L e = 30). Case E 
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FIG. 16: The potential V(R) and the energy conditions EC1= p + p r + 2pt, EC2= p + p r and 
EC3= p + pt, for 7 = 7/4, w = -3, a = 0.1 and m c = 0.1. The horizons are: r h = 7.071067814; 
r bh = 0.2003215448, r c = 4.896828668 (L e = 5); r bh = 0.2009009580, r c = 2.894500124 (L e = 3). 
Case H 
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